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ABSTRACT 

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38 /im using all modules of 
Spitzer^s IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was dis- 
covered at 17.1 /im. The feature's intensity is nearly half that of the ubiquitous 11.3 /im band. The 
7-14 /im spectral maps revealed significant variation in the 7.7 and 11.3 /im PAH features between the 
stellar ring and nucleus. Weak [Oiv] 25.9 /im line emission was found to be centrally concentrated in 
the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or 
unusually active massive star photo-ionization. Two [S III] lines fix the characteristic electron density 
in the Hii regions at Tie <200cm~'^. Three detected H2 rotational lines, tracing warm molecular gas, 
together with the observed IR continuum, are difficult to match with standard PDR models. Either 
additional PDR heating or shocks are required to simultaneously match lines and continuum. 
Subject headings: galaxies: individual (NGC 7331) - galaxies: ISM - dust, extinction - infrared: 
galaxies - techniques: spectroscopic - lines: identification 



1. INTRODUCTION 

The mid-infrared (MIR) spectra of star-forming galax- 
ies are characterized by fine structure and molecular hy- 
drogen lines, and broad emission features usually at- 
tributed to the large molecular polycyclic aromatic hy- 
drocarbon (PAH) species, set against a continuum tran- 
sitioning from the direct photospheric emission of cool 
stars at short wavelengths to the rapidly rising contri- 
bution from very small grains and heated dust toward 
longer wavelengths. The Infrared Space Observatory 
(ISO) revolutionized the use of both line and band fea- 
tures in the MIR as diagnostics of star-formation and 
AGN activity in galaxies dominated by infrared emis- 
sion, and found several key relationships between lines 
that trace excitation and the hardness of the ultravio- 
let ra diation field (e.g. [Ne ill] /[Neil], Thornle v et al 
22P0(), star-formation (e.g. 7.7 /im PAH, £gnzg^^L 
1298) and active nuclei (e.g. [Ne v] & [O iv]. lSturm et al ' 
20021^ . 



The Spitzer Space Telescope offers for the first time 
the sensitivity and spatial resolution necessary to apply 
the suite of infrared diagnostics to individual environ- 
ments in galaxies which span a full range of interstel- 
lar medium (ISM) and IR emission properties. This is 
the focus of th e Spitzer Infrare d Nearby Galaxies Sur- 
vey (SINGS — IKennicutt et alJl2?M hereafter IkoI . a 
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Legacy program for which NGC 7331 is the first galaxy 
observed by Spitzer. 

NGC 7331 has an IR luminosity typical of normal spi- 
ral galaxies with moderate star-formation activity, with 
an IR to optical ratio Lm/ Lo-pt ^ 1- By comparison, the 
extreme starburst M82 has iiR/Lopt ~ 40, and ultra- 
lumin ous infrared galaxies rang e in IR activity up to 
~400 (|Sanders & Mirabell [T99ll . Kinematic and high- 
resolution studies in both the optical and radio have been 
una ble to determine w hether the galaxy harbors an AGN 
('e.g. lFilhoetini2n?i^ . 

ISO imaging of NGC 7331 was obtained by iSmithI 
(|1998'), who found little variation in the coarse MIR spec- 
tral energy distribution traced in six broad bandpasses 
from 6-15 /im, three of which roughly match prominent 
PAH features. A surprising result from ISO was that the 
2.5-11.6/xm spectral shape is relatively invariant among 
normal galaxies, with the appealing explanation that 
the PAH molecules whose features dominate this wave- 
length region are transiently heated by single photons, 
and thus produce an emission spectrum roughly inde- 
pendent of the heating environment (Hclou ct al. 2001). 
The spectral mapping capabilities of Spitzer make it pos- 
sible to investigate variations in the PAH emission spec- 
trum within individual galactic environments, and couple 
these variations using fine structure and molecular lines 
to the physical parameters of the ISM, H ii and photo- 
dissociation regions (PDRs). 

2. OBSERVATIONS AND REDUCTION 

Spectra were obtained using all four modul es of the In- 
frared Spectrograph (IRS, Houck et al.ll2004|) in Spectral 
Mapping Mode. In this mode, the spacecraft moves in a 
rectangular raster of discrete steps, settling at each posi- 
tion before the integrations begin. For all maps, half-slit 
step sizes were used. Since Spitzer cannot roll to al- 
ter the position angle of the slits, the mapping strategy 
was designed to provide the maximum useful overlap be- 
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Fig. 1. — An overlay of the IRS footprints on the 8 fim IRAC image, with the nucleus magnified inset, highlighting the corresponding 
apertures used in LH (19-37 /^m, red) and SH (10-20 /^m, blue). The 14— 38 ^m LL radial strip is in green (atypically, following the minor 
axis due to unusual scheduling constraints for this galaxy), while the SL 5-14 /xm regions covered are in magenta and yellow (composed of 
two subslits, with the central area mapped in both subslits). The distance of NGC 7331 is 15.1 Mpc (1 kpc = 13.66"). 
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tween modules. Figure [H illustrate s these areas overl aid 
on a SINGS 8/im IRAC image fsee lRegan eraLll2?)?)l in 
this volume), and indicates the high-resolution apertures 
used for computing line parameters. 

Integration times ranged from 14-60s per pointing, 
with each position covered 2-4 times, matchi ng the uni- 
form SINGS observational strategy outlined in lK03i The 
two IRS low-resolution modules, short-low (SL, 5-14 /xm) 
and long-low (LL, 14-38 /im), provide R=50-100, while 
the two high-resolution modules, short-high (SH, 10- 
20 Aim) and long-high (LH, 19-38 /im), deliver R-600. 

The individual low-res spectra from each map w ere as- 
sembled into spectral cubes by CuBiSM f see IK03[ §6.2), 
which was also used to produce maps and extract spec- 
tra. The pipeline and calibration products used were ver- 
sion S9.5.0. Low-resolution spectra were flux-calibrated 
using the standard IRS FLUXCON tables (see Spitzer 
Observer's Manual^'^, Chap. 7), which provide unbiased 
fluxes only for point sources. We derived extended source 
flux intensities by applying three corrections: 1) an 
empirically determined aperture loss correction function 
which accounts for the pipeline's narrowing, point-source 
specific extraction aperture, 2) an estimate of the slit loss 
correction function which results from the point-spread 
function over-filling the IRS slits, and 3) an estimate 
of the integrated cross-slit beam profile. The LL back- 
ground was subtracted in situ using the extended off- 
source map wings, and the SL background was fit and 
removed using two off-source spectra. High-resolution 
spectra were extracted directly from pre-flatfield data 
and flux calibrated using stellar standards via an empir- 
ically derived relative spectral response function. This 
resulted in an absolute calibration accuracy of 25%. 

3. RESULTS 

3.1. The Spatially Varying PAH Spectrum and a 
Broad 17 ^im Feature 

The entire 5-38 /im low-resolution spectrum extracted 
from a 58"xl7" region, covering both the nucleus and a 
portion of the ring, is shown in Fig. [3 Most of the fine 
structure and II2 lines are detected in the low-resolution 
spectrum. The complete PAH spectrum from 5-20 /im is 
inset, with an untailored diffuse ISM PAH + graphite + 
silicate model of Li & Drainc (2001) overlaid. The model 
matches the four main PAH bands, at 6.2, 7.7, 8.6, and 
11.3 /im, reasonably well. The feature at 12.7 /im is con- 
taminated by [Neil] 12.8 /tm, but is also in good agree- 
ment with the model when the [Ne 11] line flux is consid- 
ered. Broad 10 /im silicate absorption is likely present, as 
the continuum between 7.7 /tm and 11.3 /im features lies 
below that inferred from bracketing band- free regions, 
most easily seen in the strong departure from the scaled 
model in this region. 

Of particular interest is the bright, broad complex at 
17 /tm, blended with H2S(1). The PAH model repro- 
duces the narrow peak on the blue wing of the feature, 
but does not predict the large width observed, which 
is approximately equal to that of the 7.7 /tm PAH band. 
Further peaks modeled at 18 and 21 /tm are not evident in 
the spectra. Deuterated PAHs are expected to have C-D 
out-of-plane bending modes which emit in the 14-16 /tm 
region, but the observed 17 /tm feature falls longward of 
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these and is much st ronger. It may be due to PAH C -C-C 
b ending modes ( e.g. lAUamandola fc Hudgin£ll2003() . 

iMoutou et alJ l)2000|) found a strong, narrow but re- 
solved emission feature at 16.4 /tm in ISO SWS spec- 
tra of the bright reflection nebula NGC 7023, the Orion 
bar, and on the surface of a Galactic molecular cloud. 
Their spectra showed no other features between 15.2 
and 17.2 /im. IRS spectra of NGC 7023 reveal a series 
of blended peaks at 16.4, 17. 4, and 17.8 ^m whic h vary 
strongly across the nebula fsee lWerner et al.l2004l in this 
volume). NGC 7331 is the first extragalactic source in 
which this feature has been discovered. 

In Fig. 131 we fit the blended 17 /tm feature as a com- 
bination of three unresolved lines matching SH identi- 
fications, and two PAH bands — one as previously re- 
ported at 16.4 /im, and a new, much stronger feature at 
17.1 /im. Removing the line contaminants, the combined 
16.4-1- 17.1 /tm feature strength is 1.5 x 10~^ Wm~^ sr~^, 
fully 49% of the 11.3 /tm feature strength, making it a 
potentially important contributor to the total MIR lumi- 
nosity, and a necessary consideration in the interpreta- 
tion of photometric redshift and luminosity distributions 
based on 24 /im MIPS or longer imaging. 

Using the MIR spectral cubes, we can test the spatial 
variability of PAH emission on kiloparsec scales. Fig- 
ure 01 illustrates three maps created from the SL spectral 
cube. The topmost map (red) is a combination of the 
7.7, 8.6, 11.3, and 12.7/tm PAH features (as indicated 
by shading on the spectra at right), and highlights the 
mapped portion of the ring and resolved nuclear core. In 
the continuum-subtracted 7.7 /tm map (green, at mid- 
dle), the nucleus, with its weak 7.7/tm PAH feature, 
nearly vanishes. In the pure continuum map at 10 /tm 
(blue, at bottom), the strong inter- feature continuum of 
the nuclear emission is demonstrated. The three spec- 
tra extracted from ring, nucleus, and both components 
together, illustrate the spatial variability of the relative 
7.7/tm and 11.3/tm PAH strengths. Although the appar- 
ent strength of the nuclear 7.7/tm PAH feature is diluted 
by excess continuum, the 7.7/tm/11.3/im ratio increases 
by ~40% from nucleus to ring. A similar band ratio in- 
crease, correlated with increasin g star formation act ivity, 
has also been observed in M31 llPagani et alJll99^ . 

3.2. H2 and Fine Structure Line Emission 

We detected seven fine structure and three molecular 
hydrogen lines, averaged over areas of size 22" x 22.4" 
(LH) and 23.1" x 15.7" (SH) (~l-1.5kpc on a side at the 
15.1Mpc distance of NGC 7331, see Fig. [T] inset). The 
fine structure lines detected arise from species ranging 
in ionization potential from 8.2 eV (Si^) to 55 eV (0^+). 
The line parameters are given in Tabled and individual 
line cut-outs are shown in Fig. |S1 

The 18.71 /im and 33.48 /tm fine structure 
lines of [Siii] are detected. For T-8000K the 
[S III] 18.71 /im/[S III] 33.48 /tm line ratio is insensitive 
to temperature and provides a density diagnostic. The 
density obtained with this diagnostic is ne<200cm~^. 
This density corresponds to a thermal pressure 
p/k < 4 X 10^ cm~'^K in the H 11 regions. 

The pure rotational lines of molecular hydrogen, espe- 
cially the S(l) and S(2) lines, which require T > 100 K for 
excitation, likely arise in PDRs surrounding or adjacent 
to these H 11 regions, and can be used to provide joint 
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Fig. 2. — The low-resolution spectrum of the inner ring and nucleus (central magenta rectangle in Fig. Q. Alternating colors indicate 

ectra from the four SL and LL orders. Inset, the expanded PAH spectrum with a diffuse ISM PAH + graphite + silicate model 

Drainc 2001) overlaid as a dashed line, arbitrarily scaled to match the 7.7 /xm peak. The broad 17 fim complex seen is blended with 

H2 S(l), and contains a sharp feature matching the model at 16.4 tiicn. F illed points represent matched photometry with filter widths shown, 

where squares denote IRAC Sum and MIPS 24 um fRegan et al.l2004l) . 
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Fig. 3. — A fit to the 17 /im feature, including unresolved [Nelll], 
[S in], and H2 S(l) in black, as well as two Drude-profile PAH com- 
ponents in blue: the previously discovered 16.4 /^m feature, and a 
new, much stronger 17.1 /xm band. 



constraints on the pressure or density and starhght in- 
tensity there. The observed H2 S(0)/S(1) and S(l)/S(2) 
hne ratios and intensities can be reproduced by PDR 
models (Draine & Bcrtoldi 2001, M.J. Kaufman et al. 
2004, in prep.) with a large fraction of PDRs in the 
beam having a similar thermal pressure as the ionized 
gas, typically T^500K and n~5000cm~'^ at the sur- 
face, with incident FUV fluxes roughly Go > 10^ times 
the local interstellar value. This is consistent with the 
strong 20-37 ^m continuum observed in Fig. 12 which in- 
dicates a component of dust with temperature 30-60 K; 
such warm dust requires Go ~ 10^-10^. 

The nucleus and inner edge of the circumnuclear ring 
in NGC 7331 therefore appear dominated by Hii regions 
and PDRs. However, our current PDR models tend to 
create too much warm dust for the observed 30 //m con- 
tinuum. This may indicate either a lack of understand- 
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'^Weighted average intensity over the entire LH map, with 4x the 
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ing of PDR heating processes (we note similar problems 
for H2 and CO rotational emission observed in dense , 
warm Galactic PDRs, see IHoUenbach &: Tielenall999D . 
or a contribution of shock heating to the H2 emission. 

Also of note is the wea k but well-detec ted [O iv] line 
near 26 ^m (see Fig. EJ. iLutz et all 1)1998^) found that 
[O iv] , with its high ionization potential of 55 eV, was 
nonetheless nearly universally present as weak emis- 
sion in galaxies classified as starbursts from the relative 
strength of their 7.7 /im PAH features. The 0''+ ion- 
ization potential is just above the He Lyman limit of 
54.4 eV, and as a result [Oiv] is rarely formed in Hii re- 
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ed by ! 
flanking regions illustrated in Fig. [Q are shown at top. In red (top panel), a sum of the main PAH features at 7.7, 8.6, 11.3 and 12.7 /xm; 
in green (second panel), the continuum-subtracted 7.7 /xm PAH feature map; and in blue (third panel) a pure continuum map centered at 
10 fim. The colored shading in the plot below illustrates the spectral regions which were used to create the three maps, including 10 /^m 
continuum subtraction in the 7.7 /im map. The spectra (bottom panel), show three separate extractions in different regions outlined with 
dashed lines in the maps: an average spectrum over the nucleus and ring (purple, lowest curve), a ring-only spectrum (orange, highest 
curve), and a nuclear spectrum (cyan, middle curve). 

The AGN hypothesis is given further credence by noting 
that prehniinary HST Paschen-a imaging of the center 
of NGC 7331 indicates httle or no atomic hydrogen Une 
emission from the nucleus (D. Calzetti et al. 2004, in 
prep.). 



4. CONCLUSIONS 

All IRS modules were used to map the inner ^^5 kpc 
of the nearby spiral galaxy NGC 7331. A new band fea- 
ture at 17.1 /xm, associated with the 16.4 /im PAH feature 
observed in Galactic nebulae, was found at 49% of the 
11.3 /im PAH strength, making it an important contrib- 
utor to the total MIR flux. Among the 7 detected fine 
structure lines was [O iv] , whose strength and central 
concentration indicate an AGN or an unusual concentra- 
tion of massive stars as the ionizing source. Two [S Hi] 
fine structure lines fix the typical pressure in H ii regions 
to be p/k<A X 10^cm~^K, and the three H2 rotational 
lines require that the bulk of the non-ionizing stellar radi- 
ation be processed through PDRs with similar pressures. 
The traditional PAH features at 7.7 and 11.3 /im were 
found to vary in strength, with the nucleus maintaining 
strong 11.3/im emission relative to 7.7/im. 
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Fig. 5. — Emission lines from the SH and LH spectra extracted 
from the inner nuclear regions (blue and red outlined) indicated in 
Fig.0 Some residual wavelength calibration offsets in the pipeline 
are evident. 



gions. Any observed [O iv] therefore requires shocks, the 
wind-enhanced far-ultraviolet emission of Wolf-Rayet 
stars, or a.n AGN, to provide su fficient high energy pho- 
tons (e.g. iSchaerer fc Stasinskai rLQOQI. 

When the [Oiv]/([Neni]-h0.44 [Neil]) ratio (0^+ to 
emissivity-weighted neon emission) is comp ared to the 
excitation measure [Ne hi] /[Neil] (as in Lut z et aP 
((T998) Fig. 2), NGC 7331 is found to lie weU above 
the locus of starburst galaxies, with [O iv] emitting more 
than 12% of the weighted neon line flux (vs. 1% typical 
for traditional starbursts). That [O iv] is centrally con- 
centrated is evident from the relative line intensity in the 
inner to outer map. Both the relative strength and tight 
spatial distribution of [O iv] support the conclusion that 
a weak AGN provides the ionizing source, though a cen- 
tral concentration of high-mass stars cannot be ruled out. 
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